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The interaction between drifting carriers and traveling electromagnetic waves is considered within 
the context of the classical Boltzmann transport equation to compute the Casimir-Lifshitz force be- 
tween media with small density of charge carriers, including dielectrics and intrinsic semiconductors. 
We expand upon our previous work [Phys. Rev. Lett. 101, 163203 (2008)] and derive in some de- 
tail the frequency-dependent reflection amplitudes in this theory and compute the corresponding 
OO \ Casimir free energy for a parallel plate configuration. We critically discuss the the issue of verifi- 

, cation of the Nernst theorem of thermodynamics in Casimir physics, and explicity show that our 

theory satisfies that theorem. Finally, we show how the theory of drifting carriers connects to 
previous computations of Casimir forces using spatial dispersion for the material boundaries. 

o ■ 

. PACS numbers: 42.50.Ct, 12.20.-m, 78.20.-i 

Q: 

CN ! I. INTRODUCTION 

• Quantum vacuum forces acting between dielectric planar surfaces or between an atom and a dielectric semi-space 
Oh| were computed long ago by Lifshitz [l[ in terms of the complex frequency-dependent dielectric permittivity e(uj) of the 

■ material boundaries. In this original formulation for ideal dielectrics, e(w) does not include contributions from current 
carriers, and as such can be called the "bare" permittivity. Extensions of the Lifshitz theory to media with large 
free charge carrier density, such as metals or highly doped semiconductors, are typically done by adding a frequency- 
q-( dependent conduction term, computed from the optical data of the material and extrapolated to low frequencies 
!— by different theoretical models (e.g., a Drude-like term iAiujq/uj, where Co is the dc Drude conductivity). The finite 
temperature Casimir-Lifshitz force is extremely sensitive to the optical response of the materials at low frequency, and 
therefore different theoretical extrapolations have resulted in conflicting conclusions about the nature of the Casimir 
force between metals and/or highly doped semiconductors. 
| Systems with small density of current carriers, such as insulators or intrinsic-semiconductors, were recently con- 

• sidered by Pitaevskii Q and by us In [IJ the thermal Lifshitz force between an atom and a conductor with 
low charge density was computed in terms of the Green function formalism, taking into account the penetration of 

, the static component of the fluctuating EM field into the conductor. This approach is quasi-static, appropriate for 
— i ■ the large distance regime of the thermal Lifshitz atom-surface interaction. The relevant (longitudinal) Green func- 
tion, expressed in terms of an auxiliary static potential field, can be computed assuming that the gas of carriers in 

• nondcgcncratc. The static potential satisfies the equation (V 2 — n 2 )(p = 0, where n 2 — 4ne 2 no/eokBT. Here — e is 
the electron charge, eo is the static bare dielectric constant of the medium (which does not take into account the 
contribution from current carriers), and no is the (uniform) carrier density. Note that k = 1/Ru is the inverse of 
the Dcbye-Hiickel screening radius i?D- For good metals the Debye radius is very small (on the order of inter-atomic 
distances), while for insulators and intrinsic semiconductors it is much larger (on the order of microns or more). 

In 0] we have extended Pitaevskii's calculation beyond the quasi-static limit and proposed a theory for the Casimir 
interaction taking into account Debye screening and carrier drift based on the classical Boltzmann equation. Rather 
than computing the force with the Green function formalism, we use the form of the Lifshitz formula written in terms 
of frequency-dependent reflection amplitudes r£ j (w) of the j-ih material boundary. Here p denotes the polarization of 
incoming waves (transverse electric TE or transverse magnetic TM) . For simplicity, we will assume that the material 
is such that there is no mixing of polarizations upon reflection (the more general case can be treated replacing the 
reflection amplitudes by 2 x 2 reflection matrices). The projection on the plane of the interface of the linear momentum 
of incoming waves is denoted by k. The Casimir-Lifshitz pressure between two plane semi-spaces separated by a gap 
of length d is 
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The prime in the sum over n means that the zero frequency n = term has to be multiplied by a 1/2 factor, and 
all reflection coefficients are evaluated at imaginary frequencies u> = i£ n , where = 2i:nk^,T/h are the Matsubara 
frequencies. 



2 



II. FROM BOLTZMANN TRANSPORT EQUATION TO REFLECTION AMPLITUDES 



In order to compute the appropriate frequency-dependent reflection coefficients for materials with small density of 
carriers we will consider that the EM field interacts with the gas of drifting carriers, and that these can be modeled as 
a continuum nondcgcnerate system. Under these conditions, it is reasonable to model the carriers with the classical 
Boltzmann transport equation coupled to Maxwell's equations for the electromagnetic field 0, Q . 

For a dielectric the carriers are charged particles (electrons or ions) hopping from site to site of the crystalline array. 
For an intrinsic semiconductor, the density of carriers and hole is equal, but their dynamics arc different; however, 
in this work we treat them as dynamically equivalent, which doubles the charge density. Assuming that there is no 
external applied field on the material, and all fields have a time dependency of the form e~ luJt , Maxwell's equations 
take the form 



V x E = j> cjH, V x H = - je(tj)cjE + J, V • E 



en 



(2) 



Here n is the carrier density, /i is the permeability of vacuum, and J = — env is the carrier current, where v is the 
mean velocity of carriers. The charge transport in the system is described by the classical Boltzmann equation, 



4- + v • V I v = -— E - — V?i - -, 

at J m n t 



(3) 



where m is the effective mass of charge carriers, vt = \fk^TJm is their mean thermal velocity, and r is the carrier 
relaxation time. Now we linearize Eq.Q with respect to the ac fields. Writing the charge density as n — > n +n(r)e _lw *, 
the current is J = — enov (after discarding the term n(r)ve~ 2lUJt ). Since by assumption there is no external field 
applied, there is no net motion of charges, Vo = 0, so the term v • Vv in Eq.([3]) behaves as e~ 2luJt , and can also be 
discarded. Finally, the linearized Boltzmann equation is 



- I v = -— E- ^Vr 

T J m no 



Inserting n = — eV • E/e into this equation, we can solve for v: 



1 — iujt 



E 



en 
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(4) 



(5) 



Plugging this expression into Maxwell's equations we derive the fundamental equation for the electric field inside the 
material, 



V- + /i e(w)w^ I + 'i 



w(l 



E 



1 + i/j e(w)- 



ujD 



ILOT 



V • (V • E). 



(0) 



Here oj c = 'iirenofj./ZJ, fi = er/m is the mobility of carriers, and D = v\t is the diffusion constant. Note that the 
frequency-dependent ratio uj c /D = 47re 2 no /l{uj)k^T coincides with k 2 = in the quasi-static limit. 

As we will show below, Eq. ([6]) allows TM and TE solutions, so that there is no cross-polarization upon reflection 
on the material, and we can safely use reflection amplitudes rather than reflection matrices in the Lifshitz formula. 
Let us assume that the material occupies the semi-space region z < and the region z > is vacuum. 



A. TM modes 



For transverse magnetic modes e y = 0, so that the electric field is 

E(r) = [e x (z)x + e z (z)z]e lkx , 



(7) 



where, from now on, we are omitting the phase factors e lu>t . Substituting this into Eq. ([6]) we obtain two coupled 
second-order differential equations for e x and e z : 



d 2 z + /ioe(w)w 1 + i— + i 



Dk 2 



3 X = ik[l + ifi i{uj)LuD}d z e z , 



-ifioe(uj)ujDd z — k + Hoe{uS)uj 1 + i — 



e z = ik[l + ifi a e(u;)uiD]d z e x , 



3 



where uj c = w c /(l — iujt) and D = D/(l — iojt). It is possible to combine these two coupled equations into two 
uncoupled fourth-order differential equations. To this end one takes the d 2 derivative of the first equation above, 
which results in terms proportional to d z e x , d 2 e x , and d z e z . This last term d z e z can be obtained from taking the d z 
derivative of the second equation above, which results in terms proportional to d 2 e x and e x . Putting all together, one 
can derive the following fourth-order differential equation for e x : 



where 



ri T (u)) = k — /ioe(Lu)uj 1 + i — 
V w 

nj (u>) = k 2 - i — 1 + i — 



(8) 



(9) 



(10) 



In a similar fashion one obtains the following equation for e z : 

(dl-ri) (d 2 z - V l)e z =0. 



A' T e r!TZ ■ 



(11) 
-ikAx/rjT- 



The solutions of these equations that vanish for z — > — oo are e x (z) = ATe nTZ + A]^e r]LZ and e z (z) 
where we assume Rc?/t > and KcrjL > 0. The amplitudes are related as A' L = —irj^AL/k and A' T 
Given this TM electric field, the associated TM magnetic field can be readily computed, H = iyATe nTZ e lkx (k 2 — 

Now we compute the TM reflection amplitude, imposing the boundary conditions on the z = interface. These 
conditions are E x , H y , and eE z continuous (the continuity of B z is automatically satisfied for TM modes). On 
the vacuum side (z > 0) the condition V • E = implies ike x + d z e z = 0, so that the fields incident on the 



interface from the z > side are E in = E [~(k z /k)± + i]e lkzZ e lkx and H in = -E ujye lkzZ e lkx / )i kc 2 , where 



we have used that in vacuum k 



W = 

arc E' 



-rE ujye' 
(A T e 



oj 2 /c 2 



The reflected fields are E r 



rE [+{k z /k)i 



ik z z e ikx an( j 



/ fiokc , where r is the reflection amplitude. The transmitted fields into the material (z < 0) 



rj T z 



A L e 
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e lkx and H* = iy{k 2 - rj^)A T e^ TZ e %kx / ii oJr} T . Im 



tit - 1 k 

posing the boundary conditions, and after some straightforward algebra, the reflection amplitude can be written as 



r(uj) = (1 — a)/(l + a), where a 
the TM reflection amplitude is 



where 



1 



k 2 - 



7c 



k 2 (k 



Expressed along imaginary frequencies lo = i£, 



1 



m)Vk 2 +e/c*-x 
mwk 2 +e/c 2 + x 

f 2 T] L r\T ~ k 21 



Vt 



(12) 



(13) 



Along imaginary frequencies, r)L and rjx take the form: 



vm - Jk 2 + 



V?nT 



(14) 
(15) 



vt{h) = ^ + e(^)g(i + m _ ( ^^ T) ) = ^k 2 + m)+^am/myc 2 , 

where er(i£) = Co/(l + £t) an d °o — e 2 noT/m are the ac and dc Drude conductivities, respectively. Therefore, Eq. 
(jT2J) gives a modified Fresnel TM coefficient due to the presence of Debye-Huckel screening and charge drift in the 
material. 



B. TE modes 



For transverse electric modes e z = 0, so that the electric field is 

E(r) = [e x (z)x + e y (z)y}e tkx . 



(16) 
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Substituting this into Eq. ((6|) we obtain two second-order differential equations for e x and 



d 2 z + fJ, e(u)u> 2 f 1 + i— + 



.bk 






- i 




e 
















e y 









0. 



e y = ik[l + i^e{uj)ujD]d z 



The solution to the first equation is e x {z) = Ae@ z , where A is a constant and /3 2 
Re/3 > 0). Plugging this solution into the second equation we obtain e v (z) = Be r,T ' 



= —ifioe(uj)ujDr]' 2 j (we assume 
+ Ce?*, where C = ikAf3{\ + 



Given this TE electric field, the associated TE magnetic field is H 

a ikx 



[l/iHouft-iBrtreir* + 



ni a e{Lo)ojD)/{p 2 

C/Je^x + Aj3eP z y + ik{Be^ TZ + Ce fiz )v\e~ 

Now we compute the TE reflection amplitude imposing the boundary conditions on the interface (the continuity 
of eE z is automatically satisfied for TE modes). On the vacuum side we have e x — 0, so that the incident fields are 
E ln = E a e lkzZ e lkx Y and H" = E (— k z x + ki)e lkzZ e lkx / ^quj, and the reflected fields are E r = rE e~ tkzZ e lkx y and 
H 1 = ri?o(+fc z x + kz)e~ lkzZ e lkx / /lout. The transmitted fields into the material are given above. A simple calculation 
leads to the expression of the reflection amplitude r = (ik z — rfr)/(ik z + rfr). Upon performing the rotation u> — > i£, 
we get 



„TE 



Vfc 2 +g 2 /c 2 -?? T 
^/k 2 +i 2 /c 2 +r] T 



;i7) 



Using Eq. (JTSJ) we see that r is the usual Fresnel TE reflection coefficient with a dielectric permittivity equal to 
the sum of the "bare" one and the ac Drude (conduction) permittivity, e(i£) = e(££) + 47rcr(i£)/£. 

As discussed in detail in 0], these modified TE and TM reflection coefficients have appropriate limiting behaviors. 
In the quasi-static limit (£ — > 0) they coincide with the ones derived in Q for conductors with small density of carriers 
in the large distance (low frequency) regime, namely rJ E (0) = (in the static limit the TE polarized field is a pure 
magnetic field, which fully penetrates the nonmagnetic material) and r™(0) = (loq — k)/(eoq + k), with q = \Jk 2 + k 2 
(in the static limit r™ interpolates between a good conductor and an ideal dielectric). On the other hand, for any 
frequency £ > 0, and in the limit of ideal dielectrics (small free charge density and small effective thermal velocity), 
we recover the usual Fresnel equations written in terms of the bare permittivity e(u). 



III. INFLUENCE OF DRIFTING CARRIERS IN THE CASIMIR-LIFSHITZ FREE ENERGY 

We now study the implications of our theory in the computation of the Casimir-Lifshitz free energy 

f = fc BT£E / ^Hl-ri,MnKMn)e- 2d ^^), (18) 

p n=0 J ^ ' 

between two identical planar semi-spaces with small density of charge carriers, such as intrinsic semiconductor media. 
As examples, we consider the cases of pure germanium and pure silicon. The reflection coefficients Eqs. (|12ll7[) that 
enter into this equation depend on temperature explicitly through the Matsubara frequencies and implicitly through 
the optical and conductivity parameters, which we proceed to quote. 

For intrinsic Ge, the bare permittivity can be approximately fitted with a Scllmcicr-typc expression 

^Heoo+^ 2 !-^, (19) 

where eo ~ 16.2, loo ~ 1-1, and usq « 5.0 x 10 15 rad/sec at T w 300K. The temperature dependence of the permittivity 
has been measured in the 20 — 300K range at wavelengths f .9 — 5.5/im Q, and shown to be very weak. Therefore, in 
this paper we assume that the permittivity is approximately constant as a function of temperature, and given by the 
above Sellmeier fitting function. The intrinsic carrier density varies with temperature as 

7i (T) = ^n v e~ , (20) 

where n c and n v are the effective density of states in the conduction and valence band, respectively. These depend on 
temperature as n c (T) = 1.98 x 10 15 T 3 / 2 cm -3 and n v (T) = 9.6 x 10 14 T 3 / 2 cm -3 (temperature is measured in degrees 
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FIG. 1: Casimir-Lifshitz free energy at T = 300K for intrinsic Ge and Si taking into account charge drift Eqs. (|12I17|I (curves 
denoted by "drift") or an additive dc conductivity term 47rcro/£ to the bare permittivity (curves denoted by "cond"). For Ge, 
(Jo = 1/(43£1 cm), and for Si, (To = 1/(2.3 x 1Q°Q. cm). The free energies are normalized to those computed using only the bare 
permittivity in the usual Lifshitz theory. 



K). The band gap energy also depends on temperature, E g (T) = 0.742-4.8 x lCr 4 T 2 /(T + 235)eV. The effective mass 
of conductivity is m = 0.12m e , where m e is the free electron mass Q. The relaxation time depends on temperature as 
t(T) = r + Tie Ci(T/300) 2 +c 2 (T/300) |j ! where Tq = .26ps, n = 1.49ps, d = -0.434, and C 2 = 1.322. At T = 300K 
one has E g = 0.66eV, n c = 1.0 x 10 19 cm~ 3 , n v = 5.0 x 10 18 cm -3 , and r = 3.9ps. 

For intrinsic Si, eo ~ 11.87, ps 1.035, and ujq ~ 6.6 x 10 15 rad/sec at T w 300K. The temperature dependence 
of the permittivity has been measured in the 20 — 300K range at wavelengths 1.1 — 5.6^tm Q, and also shown to be 
very weak. The effective density of states in the conduction and valence bands are n c (T) = 6.2 x 10 15 T 3 / 2 cm and 
n v {T) = 3.5 x 10 15 T 3 / 2 cnr 3 respectively. The band gap energy is E g (T) = 1.17 - 4.73 x 10- 4 T 2 /(T + 636)eV. The 
effective mass of conductivity is m = 0.26?7j e @- The relaxation time parameters are To = l.Ops, t\ = — 0.538ps, 
Ci = 0.0015, and C 2 = -0.09 @]. At T — 300K one has E g = 1.12cV, n c = 3.2 x 10 19 cm~ 3 , n v = 1.8 x 10 19 cm~ 3 , 
and t = 0.5ps. 

In Fig 1 we plot the Casimir-Lifshitz free energy between two identical planar intrinsic semiconducting (Gc and Si) 
semi-spaces as a function of the distance between them. We use our theory of Casimir forces with account of Debye- 
Hiickel screening and charge drift to compute the reflection coefficient Eqs. (|12I17|) and compare these predictions with 
the simple model in which the reflection coefficients are given by the usual Fresnel formulas in which the permittivity 
of the materials e(i£) is computed by adding to the bare permittivity e(z£) a dc conductivity term inao/^. In both 
models we normalize the free energies to the free energy computed using the standard Lifshitz theory using the bare 
permittivity only. For intrinsic Ge and intrinsic Si, Q c and D /£ are both very small in the relevant range of frequencies 
for the Lifshitz formula. Therefore only the n = TM mode is significantly modified by the screening and charge drift 
effects, r£ E (0) = and rJ M (0) = (e(0)t? - k)/(e(0)q + k), with q = Vk 2 + k 2 , as in Q. In all other n>l terms in Eq. 
(|18[) the reflection coefficients can be replaced by the standard Fresnel expressions in terms of the bare permittivity 
e(«£). Since for intrinsic carrier density for Ge (~ 10 13 cm -3 ) is much larger than that for Si (ss 10 10 cm -3 ), the Dcbye 
radius is of Ge, Ru = 1/n = 0.68/im is much smaller than that of Si, Rb = 24/im. As follows from Fig. 1, the 
effect of Dcbye screening and drifting carriers (denoted as "drift" in Fig. 1) becomes important for distances much 
larger than the Debye radius, so that this effect is more likely to be detected in Ge than in Si. In the latter case, for 
distances d > R^> , the Casimir force is too weak, at such a large distance, to be measured by any current or proposed 
experimental technique. From Fig. 1 we also note that when d i?rj the plates appear as perfect conductors for the 
TM n = mode, while in the case of the additive term (e(i£) + 47T(7o/£, denoted as "cond" in the figure), the plates 
appear as perfect conductors for the TM n = mode at distances of the order of At = hc/k^T (« 7^m at T = 300K), 
independent of the material properties. 
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IV. ON THE SATISFACTION OF NERNST THEOREM AS A PREREQUISITE FOR A CASIMIR 

THEORY 

Calculations of finite temperature Casimir-Lifshitz forces between media with large charge density, such as metals 
and highly doped semiconductors, have resulted in a heated debate on the adequate way to describe the optical 
properties of such systems within the Lifshitz formalism. Different phcnomenological ways to extrapolate optical data 
to low frequencies, either with a Drude model including dissipation or with a plasma model setting dissipation to zero 
from the start, result in completely different predictions for the force at finite temperature [To| . 

It has been suggested that the Nernst theorem of thermodynamics serves as a way to accept or discard conductivity 
models when applied to the computation of the Casimir-Lishitz entropy S(T) = —dE(T)/dT. The Nernst theorem 
states that the entropy of a physical system of N particles in thermal equilibrium at zero temperature is a well-defined 
constant, determined only by the degeneracy Slpf of the ground state of the system, that is, S(T = 0) = fee Infijv. 
For systems with non-degenerate ground states fijv = 1 (e.g., a perfect crystal lattice), the entropy should vanish at 
zero temperature. This is not the case for a large class of systems, including spin networks and glasses, that can have 
a large collection of degenerate ground states (with degeneracy Q,n depending on the total number of particles), so 
that S(T = 0) > in such systems. In some textbooks [ll[ it is further required as part of Nernst theorem that the 
degeneracy fijv be independent of any varying parameters of the system (such as pressure, volume, field intensities, 
etc). This has been used by some authors to discard Casimir theories which lead to a zero-temperature entropy that 
depends on the separation d between the Casimir plates. The requirement of independence of S(T = 0) on volume 
is at odds with the fact that entropy is an extensive quantity, and should grow with system size [121 ]. It is not clear 
to us that one can simply discard a model of conductivity for the Casimir plates based on the fact that the zero 
temperature entropy depends on the distance between plates. We believe that this issue requires further study. Of 
course, if for a given theoretical model for the Casimir plates S(T = 0) < 0, then such model violates Nernst theorem. 

In the remainder of this section we explicitly prove that our theory for Casimir forces with intrinsic semiconductor 
media is compatible with Nernst theorem of thermodynamics, resulting in S(T = 0) = 0, as for systems with a non- 
degenerate ground state (TIn = 1)- We will closely follow the approach in [13]. We start by expressing the Casimir 
free energy as E = A X]p, k E~=o 0g p {in8, k; 9), where 

g p (u,k;6) = ln[l-,* r£ 2 {oj, 0)e~ 2d V^ 3 ^], (21) 

and 8 = 2Trk^T /h. Note that we have allowed for an implicit dependence of the reflection coefficients on temperature. 
The Casimir-Lifshitz entropy is S = — {2ir/K)dE/dd. In Fig. 2 we plot the behaviour of g p (i£,, k; 8) as a function of the 
imaginary frequency u = i£ and as a function of k = |k| for TM and TE polarizations for different temperatures (the 
corresponding reflection amplitudes are obtained from Eqs. (| 12I17|) ) . Let us consider the TE and TM contributions 
to the entropy separately. 

For TE modes, since the reflection coefficient (|17[) depends implicitly on temperature only through Co c , which 
becomes exponentially small as low temperatures because the carrier density vanishes as T — > 0, it is possible to show 
that the 8 — > and u> — > limits of r™(w, 8) commute, that g p (u>, k; 8) is analytic in the upper-half complex lu plane, 
and that the sum over n and the derivative with respect to 8 in the expression for the entropy can be interchanged. 
Therefore, the contribution of TE modes to the entropy is [l3j 

oo' 

S TE (T) = -J2J2^ TE ( ine ^°) + m0 £ E (^8,k;8)+8gJ E ( m 8,k;8)}. (22) 

k n=0 

Here we have defined <?™ = d u g TE (u}, k; 8) and gJ M = dgg™{uj, k; 8). Using the analytical properties of the 
function g TE it is possible to write an expansion of the first two terms in (|22p in powers of temperature, re- 
sulting in (,g™(0,k)/6)r + (5^(0, k)/12)r 2 + . . ., where 5 J E (0,k) = lim^o 3<? TE (i£, k; 0)/d£ and gJ 2 E {0,k) = 
lini£_»o d 2 g(i£, k; 0)/<9£ 2 [HI]. Thus, the first two terms in (|2"2")) give a vanishing entropy at T = 0, and imply a low- 
temperature behaviour of the entropy proportional to T 2 , since <7™(0,k) = and <?™(0,k) = — e~ 2fed eg0j 2 /8fc 4 < 
(see TE plots in Fig. 2). The last term in (|22[) is proportional to dgCj c , which is exponentially small at low tempera- 
tures. Therefore, the full TE contribution to the entropy vanishes at zero temperature, namely S' TE (0) = 0. 

For TM modes, the reflection coefficient (|T2|) depends implicitly on temperature both through Cj c and D in a 
complicated fashion. Contrary to the TE case, the 8 — ► and u> — > limits of r T (u>,k;6) do not commute, and 
therefore it is not possible to write the contribution of TM modes to the entropy in the simple form (|22[) . The n > 1 
and n = terms have to be treated separately. This can be done by defining a new function g™(in8,k) which 
is identical to g™(in8, k; 8) for n > 1, and for n = it is defined as <7™(0,k) = hmg^o g™(a8,k;8), where we 
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T = 300K T = 300K 




FIG. 2: Behavior of the functions g p (i£,k) used to compute the Casimir-Lifshitz free energy and entropy for semiconductor 
materials with account of drifting carriers. The reflections coefficients are given by (|12|l and (|17[) . parameters are for intrinsic 
Ge (see text), and the distance is set to d = 1/im. The variation with temperature (in the range T — — 300K) of the TE 
function is not perceptible on the scale of the figure. The corresponding functions without account of Debye screening and 
carrier drift correspond to the T = OK plots in this figure. 



approach zero along the path uj = a9 |13[ . The TM contribution to the entropy is 

s™(T) = W .9 TM (o,M)-g TM (o,k) + £0rM (infl> k . g) + ine ~ g Tu {ine) k . e) + ^tm M] k . g)] 1 (23) 

k ^ n=0 J 

As in the TE case, the first two terms in the sum X) n =o vanish m the zero temperature limit 9^0. Their first 
non- vanishing contribution to the entropy is linear in T, since gj (0,k) > (see plots TM in Fig. 2). On the other 
hand, in the low-temperature limit the third term in the sum over n in (|23p is, 

n^o J° duJc J° dD 

The first term in pi)) is zero at T = due to the exponential decay of u c at low temperatures. Although dgD does 
not vanish at T = (since t(T) goes to a non-zero constant at T = and then -D(T) a T at low temperatures), the 
second term in (|24|) is also zero at T = because the integrand is exponentially small. Finally, in the limit 6^0, 
the first line in (|23|) vanishes since g™(0, k; 0) = g™(0, k) by definition. Therefore, the TM contribution to the zero 
temperature entropy is S™(T = 0) = 0. We conclude that our theory for Casimir-Lifshitz forces in systems with 
low density of carriers (intrinsic semiconductors, dielectrics, etc) with account of Debye screening and charge drift is 
in agreement with Nernst theorem. 



V. CONNECTION BETWEEN THIS THEORY OF DRIFTING CARRIERS AND SPATIAL 

DISPERSION 

As we have already mentioned, the static limit of the reflection coefficients (|12ll7p obtained by us using the 
Boltmann transport approach coincide with those previously derived in 0, where it was noted that the same static 
reflection amplitudes can be interpreted in terms of spatial dispersion. Indeed, using as the static permittivity 
tensor e = diag(e _L , e^, e") with the transverse permittivity e 1 - = ?o, and the longitudinal permittivity depending on 
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wavevector as e'l(fc) = €q[1 + 1/ {kRi)) 2 ], and computing the reflection coefficients for an anisotropic (uniaxial) material 
fj~3 ] , it is straightforward to recover the static versions of (jf 2|17| . 

Although the original Lifshitz paper [l[ for dispersion forces between bodies separated by vacuum did not allow for 
spatial dispersion, it can certainly be generalized to include nonlocal dielectric response in those setups, for example 
by averaging the vacuum Maxwell stress tensor and calculating field strengths via the retarded Green tensor of the 
field, that should include spatial dispersion when it is important 0, EH ■ As noted in Q , problems arise when the 
bodies are separated by a liquid instead of vacuum. 

Instead of proceeding via the Green function method, here we use the scattering formalism generalized to spatial 
dispersion to compute the force between plates separated by vacuum. The Casimir pressure between the plates is 
given by the same Eq. ([1]) , and the effects of spatial dispersion are incorporated by appropriately writing the reflection 
amplitudes r^(tu) in terms of the permittivity tensor e(w, fc) (Trl, [l7l]. In the permittivity tensor is computed in 
the random phase approximation (Lindhard dielectric function) including dissipation [l8l |. The reflection amplitudes 
are written as 



HP(k,oj) - 1 



HP(k,oj) + 1' 

where the TM and TE £f-functions are (our H, h functions correspond to the G, g functions in [l6l |) 



H™(k,Lo) 



where 



iT 1E (fc,w) 

h a (k,uj) 
h b (k,uj) 
h c (k,uj) 



-^h a (k,uj) 
h(k, to) 



( 7 o) 2 MM) 



1, 



1. 



(25) 

(26) 
(27) 



2k 



dq z 



2 7 U (M) 



2tt q 2 eW(q,uj)' 
dq z 



2tt q 2 -e^w^u/c) 2 ' 



2(u/c) 2 k-/°(k,u) 
k — 7°(fc, uj) 



dq z 



1 



-co ^ q 2 [q 2 



(28) 



Here 7°(fc,u>) = yjk 2 - {uj/cf 



= (k, q z ), and the tilde above an /i-function means the /i-function minus the same 
function except that the dielectric function is set to unity [l^]. In the limit of neglegible spatial dispersion (e" and e 1 - 
independent of q) one gets the usual Fresnel expressions for the reflection coefficients. 

We now show that our reflection coefficients can be linked to spatial dispersion even beyond the static limit, 
connecting in this way our approach based on Boltzmann transport equation with spatial nonlocality. We closely 
follow the approach of As in the static case, we assume that the nonlocal permittivity depends solely on fc, 

being independent of q z . This allows us to straightfowardly compute the integrals in (|2"8"|) . From the £f TE function 
we obtain 



H TE (k,i£) - 1 _ ^Jk 2 +i 2 /c 2 - yfk 2 + e ± (k, i£,)e/c 2 



ff TE (fc,i£) + l ^/fc 2 + e/c 2 + ^k 2 + e-L(fc, i£)e/c 2 ' 
We recover our TE reflection coefficient (JTTJ) for a transverse dielectric function independent of fc, 



-(MO = c(iO 



1 



For H™ we obtain 
where, after rotation u> 



7' 



\\(k,i 



e/c 2 



i 



k^/c 2 



k VT - Vt ^7° " (7 ) 2 



(29) 



(30) 



(31) 



Vk 2 



£ 2 /c 2 and t)t = -\/fc 

(H™ - 1)/(H™ + 1) to the expression (O for the TM reflection coefficient, one obtains H™(h,i£) 
Using this expression in (|3Tj) one can derive the longitudinal permittivity 



e x (k,i£)£ 2 /c 2 from (T^ and (J30J) . Equating 



e/c 2 



i 



ue/c 2 



km 



fc 7 ° - (7°) 2 



(32) 



As follows from the above considerations, our theory for Casimir forces with media with low density of charge carriers 
can be directly connected to spatial dispersion, as mentioned by us in 
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VI. CONCLUSIONS 

In this paper we have expanded on our previous work 0] to compute Casimir-Lifshitz forces between bodies with 
low density of charge carriers (intrinsic semiconductors, dielectrics, disordered systems, etc) taking into account 
Debyc-Hiickel screening and charge drift. Our approach is based on the classical Boltzmann transport equation, and 
is applicable to non-degenerate systems with an energy gap. We have shown how the finite conductivity of such 
systems modifies the Casimir-Lifshitz force between such materials, and made numerical predictions for the force 
using germanium and silicon plates. Our theory can be seen as a special case of spatial dispersion, and provides a 
simple way to take into account nonlocal effects in terms of readily available material properties. We have explicitly 
shown that our theory is compatible with Nernst theorem of thermodynamics. This is in agreement with previous 
work [h| that demonstrated that spatial dispersion resolves the issues with Nernst theorem. Work related to our 
approach, totally based on nonlocal dielectric responses, recently appeared [l9j ]. extending our analysis to degenerate 
systems. 

We acknowledge correspondence with C. Henkel, F. Intravaia, L.P. Pitaveskii, and discussions with S.A. Ellingsen, 
G.L. Klimchitskaya, V.M. Mostepanenko, and F.S.S. Rosa. D.A.R.D. is grateful to Victor Dodonov and the other 
organizers of the workshop on "60 Years of the Casimir Effect". He also acknowledges the support of the U.S. 
Department of Energy through the LANL /LDRD program for this work. 
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